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19.  Abstract 


This  report  summarizes  the  work  completed  on  high  frequency  seismology  by 
the  Center  for  Seismic  Studies  research  staff  between  1  October  1989  and  1 
October  1990.  The  aim  of  the  research  was  to  characterize  the  high  frequency  con¬ 
tent  of  noise  and  signals,  and  to  develop  methods  of  discriminating  between  mine 
blasts  and  other  sources  using  high  frequency  seismic  data.  Data  for  these  studies 
came  from  the  NORESS  small  aperture  array  In  Norway  and  its'  high  frequency  (up 
to  125  Hz)  central  element.  The  major  results  of  the  research  Include: 


•  The  upper  frequency  limits  of  tclcseismlc  P-wave  spectra  for  events  recorded 
at  the  NORESS  high  frequency  element  were  between  5  and  20  Hz. 

•  Path  spectra  for  explosions  from  Balapan  follow  the  attenuation  models  pro¬ 
posed  by  other  authors  (apparent  7  «0.M)  up  to  15  Hz,  but  from  15  to  20  Hz, 
the  attenuation  of  the  Balapan-NORESS  path  spectra  Is  greater  than  that  extra¬ 
polated  from  the  previous  models. 

•  Scatter  in  delay  times  between  successive  blasts  of  a  rlpple-flred  explosion 
have  a  deteriorating  effect  on  spectral  modulation. 

•  For  the  few  blasts  that  are  ideal,  i.e.,  having  no  or  negligible  errors  in  time 
delays;  the  delay  time  between  shots  as  well  as  the  number  of  shots  in  the 
explosion  can  be  estimated. 

•  The  shape  of  wavefronts  from  Balapan  explosions  as  recorded  at  NORESS  were 
nearly  identical  for  33  events  studied. 

•  Wavefronts  from  events  in  the  northeastern  section  of  the  Balapan  test  site 
were  characterized  by  a  consistent  tilt  with  respect  to  the  wavefronts  meas¬ 
ured  from  events  in  the  southwestern  section  of  the  test  site. 
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FOREWORD 


Tills  report  summarizes  the  work  completed  on  hit’ll  frequency  seismology  by 
the  Center  for  Seismic  Studies  research  staff  between  I  October  1989  and  1 
October  1990.  The  aim  of  the  research  was  to  characterize  the  high  frequency  con¬ 
tent  of  noise  and  signals,  and  to  develop  methods  of  discriminating  between  mine 
blasts  and  other  sources  using  high  frequency  seismic  data.  Data  for  these  studies 
came  from  the  NORESS  small  aperture  array  In  Norway  and  Its’  high  frequency  (up 
to  125  llz)  central  element.  A  brief  description  of  the  contents  of  the  report  arc 
given  below 

The  first  section  Is  a  study  of  the  high  frequency  content  of  teleselsmlc  P- 
waves  as  recorded  at  the  NORESS  high  frequency  element.  The  data  Included 
recordings  of  both  large  underground  nuclear  explosions  and  earthquakes.  The 
upper  frequency  limits  of  the  P-wave  spectra  for  the  33  events  studied  were 
between  5  and  20  llz.  Path  spectra  for  explosions  from  Bnlapan  follow  the  attenua¬ 
tion  models  proposed  by  other  authors  up  to  15  llz.  Prom  15  to  20  llz,  the  attenua 
lion  of  the  Balapan-NORF.SS  path  spectra  is  greater  than  that  extrapolated  from  the 

previous  models.  Other  path  spectra  follow  the  standard  attenuation  model  with 
_  • 

apparent  l  .--0.M,  but  do  not  have  data  above  15  llz. 

In  the  second  section,  data  from  the  high  frequency  element  at  NORESS  is  used 
to  estimate  the  characteristics  of  ripple-fired  explosions  in  fennoscandia.  The 
assumption  of  equal  lime  delays  between  successive  blasts  in  a  ripple-fire  leads  to 
optimistic  assesments  of  our  ability  to  evaluate  these  blasts.  The  few  case  studies 
that  have  been  undertaken  have  had  mixed  success  In  observing  any  spectral  modu¬ 
lation  effects.  We  discuss  the  deteriorating  effect  on  spectral  modulation  that 
scatter  in  delay  times  between  successive  blasts  of  a  ripple-fired  explosion  may 
have.  The  study  is  based  on  actual  observations  of  scatter  in  time  delays  of  some 
explosive  detonators  used  by  the  mining  industry.  For  the  few  blasts  that  are  ideal, 
l.c.,  having  no  or  negligible  errors  in  time  delays;  we  present  a  method  that  esti¬ 
mates  the  delay  time  between  shots  as  well  as  the  number  of  shots  In  the  explo¬ 
sion. 

The  third  and  final  section  of  this  report  is  a  study  of  slowness  and  backa- 
zimuth  estimation  using  interpolated  NORESS  data.  The  normal  25  ms  sampling 
interval  of  the  NORESS  array  elements  was  reduced  to  1  ms  through  interpolation. 
Tills  data  was  then  used  to  study  the  characteristics  of  wavefronts  from  explosions 
at  the  Soviet  test  site  Balapan,  and  provided  some  very  Interesting  results.  The 
shape  of  the  wavefronts  after  making  corrections  for  the  elevations  of  the  Indivi¬ 
dual  elements  were  nearly  identical  lor  all  events  studied.  However,  wavefronts 
from  events  in  the  northeastern  section  of  the  lest  site  were  characterized  by  a  con- 
sistant  tilt  with  respect  to  the  wavefronts  measured  from  events  In  the 
southwestern  section  of  the  test  site.  This  is  further  evidence  of  the  differences 
between  these  two  sections  of  the  test  site. 
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1 . 1  HIGH  FREQUENCY  CONTENT  OF  TELESEISMIC  P-WAVES 

ABSTRACT 

In  tills  note  wo  determine  upper  frequency  limits  of  spectra  for  P- waves  travel¬ 
ing  telescismlc  distances.  Data  recorded  at  the  111 frequency  element  of  the 
NORI'.SS  mini-array  are  analyzed  from  33  seismic  events  in  various  regions  and  of 
different  source  types.  The  procedure  to  estimate  high  frequency  spectra  using 
multiple  window  tapering  in  combination  witli  multi-sample  noise  correction 
employed  in  this  study  seems  to  extend  the  upper  frequency  limit  of  the  path  spec¬ 
tra  somewhat  compared  to  conventional  tapering  and  noise  correction  applied  to 
"low"  frequency  double  averaged  array  data.  The  P-wavc  spectra  have  upper  fre¬ 
quency  limits  between  5  and  20  II/.  'Mils  limit  Is,  to  some  extent,  dependent  on  the 
seismic  noise  that  varies  considerably  witli  daily  and  weekly  cycles.  The  highest 
upper  limits,  around  20  Hz,  were  obtained  from  some  large  explosions  in  E.  Kazakh. 
Central  Siberian  explosions  and  deep  earthquakes  in  the  Sea  of  Okhotsk  had  upper 
limits  of  about  13  Hz,  and  energy  up  to  10  Hz  were  detected  for  several  other  paths 
from  deep  earthquakes.  Spectra  shaped  by  the  path  are  constructed  by  correcting 
for  instrument  response  anti  for  source  effects.  Such  path  spectra  for  explosions  in 
the  Balapan  area  fairly  closely  follow  earlier  proposed  attenuation  models  to  the 
upper  limit  of  15  Hz.  In  the  band  15-20  Hz,  however,  the  attenuation  of  the  Balapan 
path  spectra  is  higher  than  would  be  expected  from  extrapolation  of  the  models 
beyond  15  Hz.  The  amplitude  of  path  spectra  for  other  telescismlc  paths  (  30-90°) 

tlecay  with  frequency  according  to  a  standard  attenuation  model  with  apparent 
—  * 

t  -0.14.  The  decay  is,  however,  often  over-printed  with  significant  undulations. 
Some  or  most  of  this  undulation  may  be  due  to  the  fact  that  simple  smooth  source 
models  were  used  to  construct  the  path  spectra. 
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1.1.1  Introduction 


Recordings  of  very  high-frequency  seismic  waves  have  recently  become  an 
area  of  Interest  for  discrimination  between  earthquakes  and  explosions.  A  number 
of  seismologlcal  stations  have  been  equipped  with  systems  for  recording  seismic 
waves  with  high  fidelity  up  to  and  beyond  frequencies  of  20  Hz. 

In  this  paper  we  analyze  data  collected  at  one  of  these  high-frequency  stations, 
the  NORESS  seismic  array  In  Norway,  to  study  the  propagation  of  high-frequency  P* 
waves  over  telesefsmfc  distances.  More  specifically  we  address  the  question:  at 
what  high  frequency  limit  may  significant  P-wavc  energy  be  detected  along  telese- 
ismic  paths? 

The  frequency  band  for  telesolsmic  recordings  of  P-waves  has  traditionally 
been  centered  at  frequencies  around  or  just  above  1  Hz.  However,  observations  of 
significant  energy  at  much  higher  frequencies  have  been  reported  for  some  telese- 
ismic  paths.  For  example,  Bachc  el  al.  (1985)  reported  observable  P-wave  energy 
out  to  the  high-frequency  limit  (8  Hz)  of  the  recording  system  at  the  UK  seismic 
arrays  F.KA,  CBA,  WRA,  and  YKA  from  nuclear  explosions  at  the  Eastern  Kazakh  test¬ 
ing  grounds.  The  path  from  Eastern  Kazakh  to  the  NORESS  array  has  also  been 
known  to  be  transparent  to  high  frequency  energy.  Bachc  cl  al.  (1986)  and  Walck 
(1988)  reported  significant  energy  up  to  12  Hz  and  15  Hz  for  events  at  Degelen 
Mountain  and  at  Balapan  at  the  Eastern  Kazakh  testing  grounds. 


1.1.2  Data 

T  he  data  used  here  were  recorded  at  the  high  frequency  element  of  the 
NORESS  array  in  Norway.  The  analysis  is  limited  to  data  obtained  al  a  vertical  com¬ 
ponent  Teledyne  Geolech  GS-13  seismometer  placed  in  the  surface  vault  at  the 
center  element  (site  NRAO)  of  NORF.SS.  The  data  were  sampled  at  125  samples/s 
and  the  amplitude  response  curve  of  the  instrument  and  recording  system  is 
approximately  flat  to  velocity  in  the  band  5  to  50  llz  (r/^/nre  /). 

1987  and  early  1988  data  were  retrieved  from  the  database  at  the  Center  for 
Seismic  Studies  for  the  high  frequency  element  at  NORESS.  Source  parameters  com¬ 
piled  and  reported  by  ISC  for  the  seismic  events  for  which  data  were  analyzed  are 
listed  in  Table  1.  The  paths  from  the  event  epicenters  to  the  NORESS  array  are 
shown  in  the  map  in  riflurc  2.  The  events  were  selected  on  the  basis  of  their  poten¬ 
tial  to  generate  high  frequencies  and  inc  lude  presumed  underground  nuclear  explo¬ 
sions  in  the  USSR  and  China,  deep  earthquakes,  very  large  shallow  earthquakes,  and 
earthquakes  that  were  reported  with  a  high  signal  to  noise  ratio  In  the  NORESS 
daily  bulletins.  Some  examples  of  recorded  signals  from  explosions  are  displayed 
in  sections  of  Figure  3. 
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Figure  I.  Amplitude  response  curve  (In  Volts  per  m/s)  of  the  recording  channels  of  the  high 
frequency  element  at  the  NORF.SS  array. 


TABLE  1 


EVENT  SOURCE  PARAMETERS 


Dale 

Origin 

time 

Latitude 

Longitude 

Depth 

(km) 

ttib(ISC) 

Mo(IIRVD) 

(Nm) 

1 987 

Jan 

5 

22:r>  2:40.5 

41.96 

81.32 

17 

5.9 

4.1-10'7 

1987 

Jan 

14 

11:03:48.7 

42.57 

142.85 

102 

6.5 

1 .7*1 0*9 

1987 

Mar 

3 

01:32:12.3 

46.35 

152.01 

96 

5.8 

4.6-1017 

1987 

Mar 

9 

23:58:50.4 

•25.49 

179.67 

483 

4.9 

- 

1987 

Mar 

12 

01:57:17.3 

49.93 

78.78 

1 

5.5 

— 

1987 

Apr 

17 

01:03:04.8 

49.89 

78.69 

1 

6.0 

— 

1987 

Apr 

19 

03:59:57.1 

60.5? 

57.08 

1 

4.5 

— 

1987 

Apr 

19 

04:04:5  5.G 

60.8 1 

57.55 

1 

4.4 

— 

1987 

May 

6 

04:02:05.7 

49.83 

78.13 

1 

5.6 

— 

1987 

May 

6 

04:06:14.1 

51.27 

- 1 79.90 

20 

6.3 

8.5-1018 

1987 

May 

7 

03:05:49.1 

46.74 

139.23 

430 

6.0 

1.8-1019 

1987 

May 

18 

03:07:34.1 

49.28 

147.69 

542 

6.1 

1.7-1019 

1987 

May 

19 

00:14:33.1 

29.86 

139.07 

417 

5.5 

- 

1987 

Jim 

S 

04:59:58.3 

41.58 

88.74 

1 

6.2 

2.6-1017 

198  7 

Jun 

6 

02:37:07.0 

49.86 

78.10 

1 

5.3 

— 

1987 

Jnn 

10 

14:50:1 1.5 

37.23 

21.47 

41 

5.2 

1.M017 

1987 

Jun 

20 

00:53:04.8 

49.91 

78.74 

1 

6.1 

— 

1987 

Jul 

G 

23:59:56.7 

61.50 

1 1 2.80 

1 

5.1 

- 

1987 

Jnl 

7 

18:12:53.4 

•  25.83 

178.17 

649 

5.4 

2.7*1 017 

1987 

Jul 

M 

23:46:0.3.5 

49.63 

147.83 

576 

5.9 

8.3*  1 0 17 

1987 

Aug 

2 

01:59:59.8 

73.34 

54.62 

1 

5.8 

- 

1987 

Auk 

10 

10:52:19.9 

29.87 

63.84 

165 

5.6 

1.4-1018 

1987 

Aug 

12 

01:29:56.8 

61.46 

112.76 

1 

5.0 

- 

1987 

Sep 

7 

1  1:32:27.3 

39.37 

54.76 

37 

5.5 

4.2*  1 017 

1987 

Oct 

2 

07:38:27.8 

27.35 

1 39,94 

464 

5.5 

4.9*  1 0 17 

1987 

Oct 

3 

15:14:57.4 

47.61 

56.23 

1 

5.2 

- 

1987 

Oct 

4 

18:34:22.6 

55.58 

161.62 

54 

6.0 

9.M0'7 

1987 

Oct 

6 

20:11:35.1 

52.96 

159.97 

34 

6.1 

7.31020 

1987 

Nov 

30 

19:23:19.5 

58.68 

-142.79 

10 

6.7 

7.31020 

1987 

Dec 

12 

04:51:50.5 

29.69 

140.02 

164 

6.3 

3.6’10,a 

1987 

Dec 

13 

03:21:04.8 

49.99 

78.84 

1 

6.1 

— 

1987 

Dec 

27 

03:05:04.7 

49.86 

78.76 

1 

6.1 

— 

1988 

Teh 

13 

03:05:05.9 

49.95 

78.91 

1 

0.0 

— 
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figure  2.,  The  paths  from  the  events  to  the  NOI 
here.  Asterisks,  circles  and  crosses  are  used  tc 
shallow  and  deep  earthquakes  respectively. 
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r.SS  mini-array  in  Norway  for  data  analyzed 
indicate  presumed  underground  explosions, 


Shagan  Events  at  NOR  ESS 
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Flour*  3u,  Records  from  the  NORIiSS  high  frequency  element  (vertical  component)  from  a 
presumed  underground  explosion  at  Itnlapnn  In  the  Knstcrn  Kazakh  testing  grounds.  The' 
tlmcwlndow  used  for  the  spectral  estimates  Is  Indicated  at  the  top  of  the  record  section. 
The  amplitude  scaling  differs  for  the  three  sections. 
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rigurc  3b.  Records  from  the  NORI-SS  high  frequency  element  (vertical  component)  from  a 
presumed  underground  explosion  In  Central  Siberia.  The  timewindow  used  for  the  spectral 
estimates  Is  indicated  at  the  top  of  the  record  section.  The  amplitude  scaling  differs  for  the 
three  sections. 


Ural  Events  at  NORESS 


(s) 


Floure  3c,  Records  from  the  NORI'.SS  high  frequency  clement  (vertical  component)  from  a 
presumed  underground  explosion  In  the  Ural  Mountains  region.  The  llmewlndow  used  for 
the  spectral  estimates  Is  Indicated  at  the  lop  of  the  record  section.  The  amplitude  scaling 
differs  for  the  three  sections. 


1.1.3  Data  Analysis 

An  amplitude  spectrum  was  estimated  for  the  initial  P-wave  of  eacli  recorded 
event.  The  estimates  were  based  on  a  3.0  s  time  window  to  isolate  the  first  arriv¬ 
ing  I*  wave  pulse.  The  data  window  started  0.3  s  prior  to  t he  manually  read  onset 
lime.  The  time  windows  are  marked  on  the  record  sections  in  Figure  3. 

Spectra,  calculated  from  tiie  initial  P  wave  window,  were  corrected  for  noise 
using  the  20  s  noise  segment  prior  to  the  data  window  of  the  P  signals.  In  order  to 
optimize  the  trade-off  between  leakage  resistance  and  variance  of  the  spectral  esti¬ 
mate,  we  used  the  multipie  window  spectral  analysis  method  of  Thompson  (1982). 
Tiiis  technique  has  been  applied  by  Zhu  et  at.  (1989)  to  teleseismic  recordings  of 
underground  nuclear  explosions.  After  some  experimentation  with  data,  the  time 
bandwidth  product  was  set  to  4.0  in  the  spectral  calculations. 

The  power  spectrum  of  an  initial  P  phase,  Ps (f),  was  corrected  for  noise 
according  to  the  formula: 

/\(n-/’v(/')-  PN(I). 

where  average  noise  spectra.  PN(f),  is  defined  as 


with  the  noise  spectra,  PN,( f)  (  /  l,..,(i),  being  calculated  In  the  same  manner  as  the 
spectrum  for  the  initial  P  wave  for  six  consecutive  data  windows  prior  to  the  P  wave 
window. 

Only  estimates  for  frequencies,  f,  such  that  are  used  in 

the  subsequent  analysis.  This  somewhat  conservative  cut-off  criterion  should  avoid 
any  biasing  noise  contamination. 

Single  station  spectra  are  sometimes  considered  to  be  more  ambiguous  than 
average  array  spectra  (Bachc  et  ul..  1985).  The  multiple  window  spectra  technique, 
however,  gives  smoothed  spectra  as  can  be  seen  from  the  comparison  of  spectra 
calculated  with  multiple  windows  and  with  the  10%  cosine  taper  that  is  often  used 
( Figure  A). 

The  noise  correction  applied  here  is  somewhat  different  from  those  of  Bache 
et  al.  (1985)  and  Walck  (1988),  for  example,  who  base  their  correction  on  one  sam¬ 
ple  of  noise  data  of  about  the  same  length  as  the  signal  window  that  immediately 
precedes  the  signal.  Here  we  use  an  average  of  several  consecutive  noise  samples 
(n-=6)  preceding  the  signal  and  can  thereby  reduce  the  scatter  in  the  correction  slg- 
nlficantly.  We  illustrate  this  effect  with  a  computational  exercise  on  the  noise  data 
segments.  Clearly  the  true  noise  contaminating  a  signal  power  spectrum,  Psif),  can 
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COMPARISON  OF  SPECTRAL  ESTIMATES 
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I'lQiire  •/.  Comparison  of  J’-wavc  spectra  corrected  for  instrument  response  obtained  with  a 
standard  cosine  taper  (10%)  and  with  multiple  window  tapers  for  one  of  the  events  (Apr  17, 
1087  at  Baiapan). 
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only  In;  estimated  and  not  determined  exactly.  We  can  study  the  quality  of  such  a 
noise  estimate  iisinj'  the  noise  data  segments.  Rather  than  estimating*  the  noise  por¬ 
tion  of  the  signal  data  window,  we  can  assume  that  we  want  to  estimate  the  spec¬ 
trum  of  a  noise  data  window  using  only  data  prior  to  that  window.  The  estimated 
noise  spectrum  can  then  be  compared  with  the  calculated  spectrum  for  the  actual 
noise  data  window.  Ideally  the  ratio  of  the  estimated  and  the  calculated  noise  spec¬ 
tra  should  be  equal  to  I  regardless  of  frequency.  The  standard  deviation  of  the 
ratio  can  be  used  as  a  measure  of  the  quality  of  the  noise  power  spectrum  esti¬ 
mate.  ri(iure  5  compares  the  standard  deviation  of  the  ratio  as  a  function  of  fre¬ 
quency  for  two  estimation  procedures.  The  upper  and  lower  curves  represent  esti¬ 
mates  obtained  from  one  noise  sample  with  Hanning  tapering  and  from  five  con¬ 
secutive  samples  with  multiple  window  tapering,  respectively.  The  latter  estimation 
procedure  reduces  t ho  standard  deviation  of  the  ratio  by  nearly  a  factor  of  four. 

Variations  among  spectra  based  on  recordings  at  the  same  station  from  dif¬ 
ferent  events  in  the  same  source  area  can  he  expected  due  to  differing  yields,  depth 
of  burial,  near  source  geology,  and  contribution  of  secondary  sources  like  tectonic 
release  (Bachc  et  al.,  1985).  However,  distinct  differences  among  spectra  from 
groups  of  closely  spaced  events  have  also  been  observed.  Walck  (1988)  found  that 
some  explosions  at  Sliagan  River  had  a  flat  spectral  slope  between  3  and  7  Hz, 
whereas  others  decayed  rather  smoothly  in  tills  band.  However,  signals  recorded 
under  similar  noise  conditions  from  closely  spared  sources  of  different  strength 
are  expected  to  have  different  high  frequency  limits,  and  spectral  averaging  of 
events  would  have  to  he*  based  on  weightings,  for  these  reasons  spectra  for  events 
in  Hie  same  source  region  are  not  averaged. 


1.1.4  Path  Spectra 

following  the  approach  of  Batfic  et  al.  (1985),  we  represent  the  noise 
corrected  spectra  as  path  spectra,  p(f).  The  square  root  of  the  noise  corrected 
spectrum  is  taken,  then  corrected  for  the  amplitude  response  curve,  1(f),  and  the 
source,  D(f)  : 

p(fh[h(f)f/(i(f)-mn). 


Different  models  have  been  used  in  the  literature  for  source  corrections,  D(f). 
Bachc  et  al.  (1980)  and  Walc  k  (1988)  use  a  model  by  Mueller  and  Murphy  (1971), 
and  /.hit  et  al.  (1989)  use  a  model  by  Von  Seggern  and  Blnndford  (1972).  Both  of 
these  models  are  mainly  derived  from  observations  al  the  Nevada  test  site.  Their 
validity  for  the  events  studied  here  with  regard  to,  for  example,  near  source  struc¬ 
ture  is  unknown.  Furthermore,  there  was  no  information  available  on  parameters 
like  shot  depth  and  yield  for  t he  explosions  in  Table  1.  Because  of  these  uncertain¬ 
ties  we  employ  the  somewhat  simpler  model  of  Von  Seggern  and  Blnndford  (1972) 
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lifjun’  r>.  Standard  deviation  of  the  ratio  of  estimated  and  at tit  a  I  noise  power  spectra.  The 
estimated  spectrum  for  the  upper  curve  is  based  on  one  data  window  immediately  preceding 
the  data  window  for  the  actual  noise.  Standard  Hanning  tapering  was  applied  in  the  spectral 
rah  illations.  The  estimated  spectrum  for  the  lower  curve  represents  an  average  of  spectra 
lor  five  data  windows  immediately  preceding  the  actual  noise.  Multiple  window  tapering 
was  applied  In  the  calculations.  The  two  curves  illustrate  that  the  combined  use  of  several 
noise  data  windows  and  multiple  window  tapering  reduces  the  standard  deviation  In  the 
estimate  of  the  actual  noise  significantly. 


0  10  20  30  40  50 

FREQUENCY  (HZ) 


M2 


for  a  granitic  source  medium: 


[A-'wr/kVi  if 

[(2*r/k)’+  l)V' 


The  parameters  /1=4.08  and  k~28.727frW'  The  yield,  W,  is  assumed  to  scale  with 
m/,(/.9C)  according  to  the  linear  relation  mj,(/.VC)=4.55  1 0.7log)0(W)  estimated  by  Ver- 
gino  (1089).  The  iny  -  yield  relation  is  based  on  announced  yield  Information  for 
explosions  in  Eastern  Kazakh.  Reduced  displacement  potentials  calculated  for  this 
model  are  shown  in  Figure  6  for  the  magnitude  range  O=*4.5-6.0. 

Due  to  the  relation  between  magnitude  and  corner  frequency  at  high  frequen¬ 
cies,  the  source  term  can  be  simplified  for  events  above  a  threshold  magnitude. 
The  body  magnitudes,  mi, (ISC’)  of  the  explosions  in  Table  1,  are  all  equal  to  or 
greater  than  4.5.  Bache  el  al.  (1985)  concluded  that  source  corner  frequency  effects 
are  contaminating  path  spectra  for  K.  Kazakh  explosions  below  2.5  Hz  or  so  for 
events  with  mh(lSQ>S.A.  Bache  el  al.  (1986)  find  that  source  corner  effects  can  be 
seen  up  to  3  Hz  for  the  small  explosions  at  Degelen  Mountain  (  my(ISC)- 4.5),  The 
curves  in  Figure  6  for  the  model  by  Von  Seggern  and  Blandford  (1972)  show  source 
coiner  effects  around  5  llz  for  my(ISQ- 4.5.  This  suggests  that  the  path  spectra 
obtained  with  Dif^/f2  are  valid  at  frequencies  greater  than  about  5  Hz. 

lor  earthquakes,  we  use  l)(f)~l/f*.  Table  I  includes  estimates  of  seismic 
moments,  M„,  of  the  earthquakes  (IIVKD  solutions  published  in  the  ISC  Monthly  bul¬ 
letins).  We  assume  a  standard  omega  square  model  for  the  earthquake  sources,  i,e„ 
M„/(  I  i  (f/fc)7),  with  the  corner  frequency  /'.  scaling  with  seismic  moment  as 
Using  corner  frequency  and  moment  estimates  for  some  50  earth¬ 
quakes  in  California,  Turkey,  and  Tonga-Kcrmndet  (Hanks  and  Boore,  1984;  Molnar 
and  Wyss,  1972),  we  normalize  the  constant  fc«(1.77*I0,<),^J  for  fc~  1.0  Hz.  Accord¬ 
ing  to  this  corner-frequency  model,  the  smallest  seismic  moment  of  the  earthquakes 
studied  would  correspond  to  a  corner  frequency  of  0.5  llz.  This  means  that  a 
source  correction  of  the  form  D(f)-\/(?  would  be  sufficient  for  calculating  path 
spectra  for  these  earthquakes. 


Hal  a  pan 


Figure  7  shows  the  path  spectra  for  six  Balapan  events  that  are  compared  with 
the  hybrid  model  estimated  by  Walck  (1988)  for  Balapan.  This  model  assumes  that 
the  attenuation  as  a  function  of  frequency.  C(f),  can  be  written  as: 


C(f)~ex  p 


2*  /*•/,'> -tan  *• 


I  |  4  S  f2TMTm 
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Von  Seggern-Blandford 
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ligurtt  6.  lar  field  displacements  as  a  function  of  frequency  for  an  explosion  source  model 
by  Von  Severn  and  Blandford  (11)72)  for  m/,  between  4.5-G.0.  The  curves  are  normalized  to 
1  at  zero  frequency  to  compare  the  relative  location  of  overshot  and  corner  frequency  as  a 
function  of  source  strength. 
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SHAGAN  RIVER 


['lyure  7.  Path  spectra  for  llalapan  test  site  explosions.  An  attenuation  model  by  Walck 
(1988)  is  drawn  through  the  spectra  at  5  II/.,  and  a  standard  attenuation  model  T  »0.14  Is 
drawn  in  the  lower  part  of  the  figure.  Hie  amplitude  scale  is  arbitrary  since  the  spectra  have 
been  shifted  for  clear  display. 
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with  parameter  values  f<‘,-0.6  s,  rm-0.05  s,  t\MU2  s,  and  rM« 10000  s.  Three  of  the 
events  (87/03/12.  87/04/17,  and  87/06/20)  were  analyzed  by  Walck  (1988)  for  the 
NOR  ESS  array  using  data  sampled  at  40  samples/s. 

All  path  spectra  follow  G(()  fairly  closely,  l.e.,  p(/‘)«C(/‘)»  lip  to  some  fre¬ 
quency  between  13  and  15  Hz.  However,  beyond  13-15  Hz  they  drop  off  faster  than 
G(f),  i.e.,  p(f)<G(f).  Clear  differences  among  the  path  spectra  can  also  be  noted. 
The  data  for  two  of  the  events  (87/12/13  and  87/12/27)  closely  follow  the  hybrid 
model  up  to  15  Hz,  above  which  the  attenuation  Increases  significantly  compared  to 
the  extrapolated  curve  for  the  hybrid  model.  Both  spectra  have  significant  energy 
up  to  and  slightly  above  20  Hz.  The  data  for  two  other  events  (87/04/17  and 
88/02/13)  also  go  up  to  about  20  Hz.  These  path  spectra  seem  to  decay  a  bit  faster 

than  the  hybrid  model,  hut  thcy.  like  the  other  pair  of  path  spectra,  also 

drop  in  amplitude  at  around  15  Hz.  The  remaining  two  events  (87/03/12  and 
87/06/20)  contain  significant  energy  up  to  about  15  Hz;  and  they  follow  the  hybrid 
model  up  to  around  13  Hz,  above  which  they  drop  below  the  model. 

The  relative  locations  of  the  Balapan  events. are  shown  in  figure  8.  There  is  no 
obvious  correlation  between  shapes  of  the  patli  spectra  and  location. 


Degelen  and  Central  Siberia 

Figure  9  shows  spectra  for  two  Degelen  events  and  three  Central  Siberian 
events.  These  events  were  also  analyzed  by  Walck  (1988)  and  Walck  (1989)  using 
the  NORF.SS  40  Hz  array  data. 

The  estimated  path  spectra  are  compared  with  the  hybrid  model  for  Degelen 
estimated  by  Walck  (1988)  (with  parameter  f,*,-0.6  s,  rw=0.05  s,  fi=0.12  s,  and 
Tsi-  10000  s).  The  straight  line  corresponds  to  a  simple  attenuation  model, 
G(f)^cxp(-irft‘)  with  slope  t‘^0.14,  representing  Ihc  apparent  t*  value  (Der  et  a!., 
1986). 

The  data  for  the  two  Degelen  events  are  in  reasonable  agreement  with  Walck’s 
results  and  her  estimated  hybrid  model.  However,  ihc  spectra  here  go  up  to  15  Hz 
compared  to  13  Hz  for  Walck  (1988). 

Similarly  the  spectra  for  the  the  Central  Siberian  events  go  to  somewhat  higher 
frequencies  (13  compared  to  12  11/  for  Walck,  1989).  furthermore,  they  seem  to 
decay  somewhat  less  than  predicted  by  a  t‘~0. 14  attenuation  obtained  by  Walck 
(1989).  This  may  be  due  to  the  fact  that  different  source  models  were  used  In  cal¬ 
culating  the  path  spectra.  The  events  are  rather  small  (  mb(ISQ» 5)  and  there  may 
be  some  source  spectrum  corner  frequency  effec  t  up  to  about  5  llz. 
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figure  9.  Path  spectra  for  two  explosions  In  Mount  Dcgelcn  and  three  closely  spaced  explo¬ 
sions  in  Central  Siberia.  An  attenuation  model  by  Walck  (1988)  Is  drawn  through  the  spectra 

—  * 

at  5  Hz,  and  a  standard  attenuation  model  t  ^O.H  is  drawn  In  the  lower  part  of  the  figure. 
The  amplitude  scale  is  arbitrary  since  llu*  spectra  have  been  shifted  for  clear  display. 
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Slnkiang-Ka/.akh-Urals-NOliFSS  Traverse 


The  paths  for  events  at  the  testing  {’rounds  In  S.  Sinkiang,  E.  Kazakh,  and  the 
two  events  in  the  Ural  Mountains  are  almost  at  the  same  azimuth  from  NORF.SS  (Fig- 
lire  2.).  The  path  spectra  for  these  explosions  arc  compared  In  Flyura  10. 

The  hybrid  model  by  Walck  (1088)  for  lialapan  Is  drawn  for  comparison 
together  with  a  standard  attenuation  model  with  /‘r-0.14. 

The  spectra  of  the  two  Ural  Mountains  events  are  very  similar  and  follow  the 
Walck  hybrid  model  fairly  closely  in  the  band  3-1 4  llz,  although  they  appear  some¬ 
what  flatter  than  the  model  at  frequencies  below  10  llz.  It  should  be  noted  that  the 
Ural  Mountain  events  are  much  smaller  in  magnitude  and  the  source  corrections 
applied  here  may  not  be  appropriate.  The  path  spectrum  for  S.  Sinkiang  contains 
much  less  high  frequencies  than  the  Walck  model.  In  fact,  It  seems  to  decay  even 
faster  than  an  apparent  f*«0.M.  Flyura  II  shows  the  raypaths  in  the  crust  and  the 
upper  mantle  from  these  sites  to  NOKESS.  The  paths  from  Ural  Mountains  and  E. 
Kazakh  bottom  just  above  the  600  and  1500  km  discontinuities,  respectively;  and 
the  path  from  Sinkiang  goes  just  underneath  the  1500  km  discontinuity.  A  Q  model 
by  Der  et  al.  (1986)  has  a  low  Q  layer  between  100  and  200  km  (LVZ,  in  Figure  11) 
and  a  sharp  increase  of  Q  at  600  km  depth.  There  is  no  obvious  correlation 
between  this  Q  model  and  the  fact  that  the  path  spectra  for  the  Sinkiang  event  Is 
depleted  in  high  frequency  energy  relative  to  those  of  the  Ural  Mountains  and  H. 
Kazakh  events. 


Deep  Farthqaakes 

Figure  12  shows  path  spectra  for  the  deep  earthquakes  analyzed  here.  Most  of 
them  contain  significant  energy;  up  to  10  llz.  Two  events  In  the  Se.  of  Okhotsk 
(w/,(/5C)« 6.0)  contain  frequencies  up  to  about  13  llz.  Most  of  the  data  follow  a 
/*=0.M  trend;  however,  but  there  Is  clear  variations  around  this  trend  in  the  path 
spectra. 


Shallow  Earthquakes 

The  three  shallow  earthquakes  (m(ISQi,~ 5.5-6.2)  in  Figure  13  contain  signifi¬ 
cant  frequencies  up  to  about  8  llz,  but  the  spectra  show  quite  a  bit  of  variation 
around  the  attenuation  (*-0.M. 
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rijurt’  in.  Path  spectra  for  the  traverse:  Slnklans'-r.Kazakh-IJral  Mountains-  NORI-SS.  An 
attenuation  model  by  Walck  (l?)88)  Is  drawn  Ihrnunh  the  spectra  at  5  II/.,  and  a  standard 

_ a 

attenuation  mode!  t  =0.M  Is  drawn  in  the  lower  part  of  the  figure.  The  amplitude  scale  Is 
arbitrary  since  the  spectra  have  been  shifted  for  clear  display. 
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Figure  11.  Schematic  of  ray  paths  for  explosions  along  the  traverse:  Sinkiang-K.Ka7akh-Ural 
Mountains-  NORESS.  Major  discontinuities  are  indicated  at  100-200  (marked  I.VZ  and  low  Q 
7one  according  to  Dor  el  al„  1080).  000  and  1500  km. 
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Deep  Focus  Events 
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ritjurc  12.  Path  spectra  for  deep  focus  earthquakes.  A  standard  attenuation  mode)  f  =0.14 
is  drawn  through  the  spectra  at  S  II/..  The  amplitude  scale  is  arbitrary  since  the  spectra  have 
been  shifted  for  clear  display. 
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Fiyure  lJ.  Path  spectra  of  some  shallow  earthquakes.  A  standard  attenuation  model 
1  =0.14  is  drawn  through  the  spectra  at  5  Hz.  The  amplitude  scale  is  arbitrary  since  the 
spectra  have  been  shifted  for  clear  display. 
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Regional  Distances 


l  idiire  14  shows  path  spectra  for  events  at  distances  less  than  30®.  The  explo¬ 
sion  spectra  have  been  corrected  for  source  by  the  Von  Seggern  and  Blandford 
( I  iJ 72)  model,  and  the  earthquake  spectrum  has  been  multiplied  by  a  factor  of  /*. 

The  events  in  the  Ural  Mountains  show  much  less  attenuation  than  the  other 
events,  which  approximately  follow  a  simple  f’--0.14  model. 


PKP  phases 

Figure  15  shows  the  path  spectra  for  two  deep  focus  events  south  of  the  Fiji 
Islands  (  m(lSC)i, =4.9  and  5.4)  at  an  epicentra!  distance  of  143.9®.  There  Is  signifi¬ 
cant  energy,  up  to  almost  10  II/,  for  the  larger  of  the  two  events. 


1.1.5  Discussion 

The  path  spectra  for  the  F..  Kazakh  events  are  strikingly  different  from  those  of 
other  events  at  teleseismic  distances.  They  have  the  highest  upper  frequency  limits, 
between  15-20  Hz. 

The  upper  limit  of  the  amplitude  spectrum  is,  of  course,  to  some  extent  depen¬ 
dent  on  the  prevailing  seismic  noise  level,  which  varies  significantly  in  a  cyclic 
manner  with  time.  This  is  illustrated  by  Figure  16.  which  shows  the  range  of  varia¬ 
tion  of  noise  power  for  the  data  studied  here.  E.  Kazakh  explosions  are  usually  set 
off  early  in  the  morning,  local  NORESS  time,  when  the  noise  level  is  low. 

The  E.  Kazakh  spectra,  however,  appear  exceptional  not  only  with  regard  to 
upper  frequency  limits,  but  more  importantly  with  regard  to  shape.  They  have 
much  less  amplitude  decay  with  frequency  than  those  of  other  events  at  teleseismic 
distances. 

The  path  spectra  of  the  two  Ural  Mountain  explosions  at  regional  distance  and 
at  virtually  the  same  backazimuth  as  the  E.  Kazakh  explosions  are  quite  similar  to 
the  E.  Kazakh  spectra.  For  the  Sinkiang  explosion,  also  at  the  same  backazimuth  but 
with  path  lengths  greater  than  the  E.  Kazakh  events,  an  apparent  significant  deple¬ 
tion  of  high  frequency  energy  in  the  path  spectrum  is  observed.  One  might  specu¬ 
late  that  a  high  Q  zone  extends  from  NORESS  to,  and  perhaps  somewhat  beyond,  E. 
Kazakh.  It  should  be  noted,  however,  that  the  Ural  Mountain  path  spectra  are  highly 
dependent  on  the  assumptions  of  the  source  spectra  corrections  because  of  the 
small  explosion  sizes.  The  sensitivity  of  path  spectra  to  assumptions  "about  source 
spectra  is  illustrated  by  results  for  the  three  Central  Siberian  explosions.  The  path 
spectra  by  Walck  (1989)  using  the  source  functions  of  Mueller  and  Murphy  (1971) 
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Figure  14.  Path  spectra  for  some  events  at  regional  distances  (i.e.,  less  than  30°).  A  stan- 
_  * 

dard  attenuation  model  t  -0.14  is  drawn  through  the  spectra  at  5  Hz.  The  amplitude  scale 
is  arbitrary  since  the  spectra  have  been  shifted  for  clear  display. 
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I  iyurc  15.  Path  spectra  for  PKP  phases.  A  standard  at  tent  nation  -iodel  t  =0.14  is  drawn 
thtongh  the  spectra  at  5  Hz.  The  amplitude  scale  is  arbitray  since  the  spectra  have  been 
shifted  for  clear  display. 
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Figure  16.  Range  of  variation  for  noise  power  spectra  for  the  data  studied  here.  The  lower 
and  upper  curves  represent  minimum  and  maximum  power  taken  over  all  individual  spectra. 
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closely  follow  the  standard  model,  f  -0,14,  whereas  the  spectra  obtained  here  with 
the  source  function  of  Von  Seggern  and  fllandford  (1972)  appear  to  be  attenuated 
less  than  the  standard  model. 

The  anomalous  spectral  enrichment  of  I7..  Kazakh  path  spectra  have  been 
pointed  out  earlier.  Ilache  ct  al.  (1985)  found  that  the  high  frequency  attenuation 
in  recordings  at  NORSAR  for  I7.  Kazakh  events  were  quite  different  from  that  in 
recordings  at  UK  arrays.  They  suggested  that  characteristics  of  the  near  receiver 
portion  of  the  path  to  NORSAR  must  be  important  and  that  the  difference  for  the 
NORSAR  path  appeared  to  be  associated  with  some  difference  In  the  contribution  of 
scattering  to  attenuation  of  the  initial  P-wave. 

Walck  (1989)  argues,  however,  on  the  basis  of  the  differences  In 
NORF.SS/NORSAR  and  in  UK  recordings  of  I'.  Kazakh  events  that  a  region  of  low 
attenuation  can  not  lie  located  near  the  source  region  and  that  one  could  not 
further  constrain  its  location  along  the  remainder  of  the  path  between  Eastern 
Kazakh  and  NORESS. 

_* 

Walck  (1989)  goes  on  to  conclude  from  an  estimated  'normal'  (i.e.,  t  =0.14) 
shield  attenuation  for  the  path  from  Central  Siberia  to  NORESS  that  since  change  in 
source  region  produces  dramatically  different  path  spectra  at  NORESS,  the  receiver 
structure  immediately  beneath  NORESS  and  NORSAR  is  probably  not  responsible  for 
the  anomalous  attenuation  behavior  for  the  E.  Kazakh  -  NORESS  path,  and  that  a 
high  Q  zone  must  be  located  elsewhere  along  the  path. 

The  study  of  Q  is  in  its  infancy,  and  results  on  its  variation  with  frequency  and 

its  radial  and  lateral  variations  are  afflicted  with  considerable  uncertainties  (Ander- 

_  * 

son,  1989).  Der  el  al.  (198f>)  obtain  t  values  from  a  tentative  Q-model  for  the 
shield  areas  of  Eurasia  at  a  distance  of  60*  that  vary  between  0.11-0.21  In  the  fre- 
quency  band  1-10  llz.  Several  of  the  path  spectra  decay  largely  like  t  =0.14  regard¬ 
less  of  epicentral  distance,  at  least  in  the  teleseismic  window.  This  is  in  qualitative 
agreement  with  the  results  of  Der  et  al.  (1986)  and  with  Q  increasing  with  depth. 

Even  if  the  path  spectra  for  several  events  of  different  source  types  decay 
largely  like  l  =0.14,  the  decay  is  often  over  written  by  significant  undulations. 
Some  or  most  of  this  undulation  may  be  due  to  the  fact  that  the  simplistic  source 
models  used  to  construct  the  path  spectra  do  not  adequately  represent  the  actual 
source  spectra.  Ibis  means  that  the  calculated  patli  spectra  are  over-printed  by 
source  effects  that  may  not  readily  be  removed.  This  may  make  it  difficult  to  apply 
source  models  that  are  smooth  as  a  function  of  frequency  to  estimate  source  scal¬ 
ing  laws. 
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1.1.6  Concluding  Comments 

In  tills  note  wc  determine  upper  frequency  limits  of  spectra  for  P-waves  travel¬ 
ing  over  teleseismic  paths.  Data  recorded  at  the  high  frequency  element  of  the 
NORESS  mini-array  are  analyzed  from  33  seismic  events  In  various  regions  and  of 
different  source  types. 

The  high  frequency  data  and  the  procedure  to  estimate  spectra  using  multiple 
window  tapering  in  combination  with  multi  sample  noise  corrections  employed  here 
seem  to  extend  the  upper  frequency  limit  of  the  path  spectra  somewhat  compared 
to  conventional  tapering  and  noise  corrections  applied  to  "low"  frequency  array 
data. 


The  P-wave  spectra  have  upper  frequency  limits  between  5  and  20  Hz.  This 
limit  is.  to  some  extent,  dependent  on  the  seismic  noise  that  varies  considerably 
with  daily  and  weekly  cycles.  The  highest  upper  limits,  around  20  Hz,  were 
obtained  from  some  large  explosions  in  E.  Kazakh.  Central  Siberian  explosions  and 
deep  earthquakes  in  the  Sea  of  Okhotsk  had  upper  limits  of  about  13  Hz,  and 
energy  up  to  10  Hz  was  detected  for  several  other  paths  from  deep  earthquakes. 

Spectra  shaped  by  the  path  are  constructed  by  correcting  for  instrument 
response  and  for  source  effects.  Such  path  spectra  for  explosions  at  the  Balapan 
test  site  area  follow  earlier  proposed  attenuation  models  fairly  closely  up  to  the 
upper  limit  of  15  Hz.  In  the  15-20  Hz  band,  however,  the  attenuation  of  the  Balapan 
path  spectra  is  higher  than  would  be  expected  from  extrapolation  of  the  models 
beyond  15  Hz.  This  illustrates  that  extrapolation  outside  the  frequency  band  of 
actual  recording  may  give  misleading  results.  Actual  observations  at  high  frequen¬ 
cies  are  important  in  assessing  the  value  of  high  frequency  data  for  seismic  verifi¬ 
cation. 

Path  spectra  for  other  teleseismic  paths  (  30  90*’)  are  reduced  as  a  function  of 
frequency  according  to  a  standard  attenuation  model  with  apparent  t  =0.1-4. 

The  reduction  of  these  path  spectra  is,  however,  often  over-written  by  a  signi¬ 
ficant  undulation.  Some  or  most  of  this  undulation  may  be  due  to  the  fact  that  sim¬ 
ple  smooth  source  models  were  used  to  construct  the  path  spectra.  This  undulation 
may  not,  therefore,  be  easily  removed. 


Hans  Israelsson 
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2.1  ESTIMATING  CHARACTERISTICS  OF  RIPPLE-FIRED  EXPLOSIONS 


2.1.1  Introduction 

Identifying  large  numbers  of  chemical  explosions  for  mining  and  other  indus¬ 
trial  purposes  presents  a  major  challenge  to  seismic  verification  systems  operating 
in  the  low  magnitude  range:  2.0-2.5.  In  a  study  of  mine  blast  activity,  Richards  et  al. 
(1990)  estimated  that  in  1987  there  were  about  15,000  chemical  explosions  in  the 
U.S.,  virtually  all  ripple-fired,  with  charges  larger  than  50  tons.  The  data  in  figure  1, 
showing  ML  magnitudes  and  explosion  charge  for  some  ripple-fired  explosions  in 
Scandinavia,  indicate  that  magnitudes  between  2-2.5,  indeed,  correspond  to  charges 
of  50-100  tons. 

Some  methods  of  identifying  ripple-fired  explosions  draw  upon  the  assump¬ 
tion  of  equal  time  delays  between  successive  blasts  in  a  ripple-fire.  In  theory,  spec¬ 
tra  of  recorded  signals  will  then  be  overprinted  by  a  modulation  pattern  that 
should  be  easily  observed  and  that  directly  relates  to  the  time  delay  between  blasts 
(see  e.g.,  Baumgardt,  1988).  Case  studies  with  these  methods  have,  so  far,  been  lim¬ 
ited  to  small  numbers  of  events  compared  to  the  scale  of  the  overall  activity  of 
chemical  explosions  (Baumgardt,  1988;  Smith,  1989;  Hedlin  et  al.,  1989;  Kim  et  al, 
1990).  Furthermore,  these  studies  appear  to  have  had  mixed  success  in  observing 
any  spectral  modulation  effects.  In  some  cases,  accurate  and  independent  informa¬ 
tion  on  the  characteristics  of  tiie  mining  explosions  is  not  available  to  confirm  con¬ 
clusions  from  observed  spectral  modulation. 

In  this  note  we  discuss  the  deteriorating  effect  on  spectral  modulation  that 
scatter  in  delay  times  between  successive  blasts  of  a  ripple-fired  explosion  may 
have.  The  study  is  based  on  actual  observations  of  scatter  in  time  delays  of  some 
explosive  detonators  used  by  the  mining  industry.  In  addition  we  present  a  method 
attempting  to  estimating  the  autocorrelation  for  the  modulation  function  of  an  ideal 
ripple-fired  explosion,  i.e.,  having  no  or  negligible  errors  in  time  delays.  This  some¬ 
what  novel  approach  allows  an  estimate  of  not  only  the  delay  time  between  shots, 
bu»  also  the  number  of  shots  of  the  explosion. 


2.1.2  Deterioration  Of  Modulation  Function  With  Scatter  in  Delay  Times 

The  standard  model  for  the  source  function  of  a  ripple-fired  explosion,  x(t), 
consisting  of  N  equal  explosions  delayed  by  the  inlet  vai,  t,  assumes  that  it  can  be 
written  as: 

N-\ 

X(l)~  i,'  S(t-k-r). 

K() 
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liiline  1.  Local  magnitudes,  Ml,  plotted  against  charges  of  ripple-fired  explosions 
for  two  paths  in  I'ennoscandiu.  Information  on  explosions  have  been  obtained  from 
a  database  compiled  by  Dahle  el  al.  (1989).  Crosses  represent  ML  in  the  NORLSS 
bulletin  from  dynamex  explosions  in  the  Dala  Kalk  open  pit  (limestone).  Open  cir¬ 
cles  aie  for  ML  al  the  station  SO  of  the  Finnish  Station  Network  from  slurry  explo¬ 
sions  in  the  Sydvaranger  open  pit  mine  (iron  ore).  The  nominal  charges  of  slurry 
explosions  have  been  multiplied  by  a  grade  factor  of  0.9. 
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Here  5(f)  is  the  source  function  for  each  individual  explosion.  By  applying  the 
Fourier  transform  to  the  source  function  for  the  total  ripple-fired  explosion,  the 
source  function  for  an  individual  explosion  and  a  modulation  function  are 
separated: 

X(jw)~S(jw)-M(j<.o). 


The  Fourier  transform  of  the  modulation  function,  can  in  turn  be  represented 

as: 

N--I 

i:  expi-j-wk  r). 

k  rtO 


If  we  assume  that  there  is  a  scatter  in  the  delay  time.  rk,  for  shot  k,  we  can  write: 
Tk=k’T-\-ik,  where  the  error  term  ik  is  assumed  to  be  normally  distributed  with  zero 
mean  and  standard  deviation  ak.  The  detonators  of  a  ripple-fired  explosion  go  off 
independently  and  tk  are  therefore  independent  (Fauske,  1990).  Moreover,  it  is  rea¬ 
sonable  to  assume  that  the  standard  deviation  increases  linearly  with  shot  number 
(rauske,  1990  and  Stump  and  Reamer,  1990),  e.g.,  nk=o()+kir\.  The  average  value 
(expectation  operator)  of  each  exponential  term  in  the  sum  that  makes  up  the 
Fourier  transform  of  the  modulation  function,  M(/w)  is: 

i:[exp(-j-uik-T\(k))]==i'X[)(-j-u-rk-  J^-w2-(<r„  hA'fr,)2). 


It  is  not  obvious  how  to  reduce  the  expression  for  the  average  of  the  power  spec¬ 
trum  of  the  modulation  function  into  a  simple  analytical  form  because  k,  the  shot 
number,  is  quadratic.  However,  the  average  of  the  power  spectrum  of  the  modula¬ 
tion  function,  M(/w)-M*(/w),  can  be  written  in  a  simple  closed  form  if  <T|=0  (Israels- 
son,  1990). 


Data  on  Scatter  in  Delays 

Estimates  of  the  standard  deviation,  ok,  as  a  function  of  shot  number  have 
been  given  for  some  detonators  used  on  the  Norwegian  market  by  rauske  (1990). 
Examples  of  these  estimates  are  shown  in  Figure  2,  where  the  standard  deviation  is 
plotted  as  a  function  of  percent  of  the  delay  time,  r,  for  millisecond  (  7-25  ms)  and 
quarter  second  (r^250m5)  detonators.  For  example,  this  means  that  for  the  25  ms 
delay  times  the  standard  deviation  increases  from  about  2  to  6  ms  for  shot  numbers 
k-~  1  and  k~ 20,  respectively. 

These  detonators  are  used  for  no  more  than  25  intervals,  i.e.,  k<25,  and  the 
standard  deviations  are  clearly  within  the  requirement,  crk<r/3,  that  such  detonators 
are  supposed  to  meet  (Fauske,  1990).  This  requirement  would  avoid  successive 
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STANDARD  DEVATfON  (par  ce 


STANDARD  DEVIATION  IN  DELAY  TIME 


Figure  2.  Standard  deviation  of  time  delay  (in  percent  of  delay  between  successive 
shots)  as  a  function  of  shot  number  (al  ter  I'auske,  1990). 
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shots  being  fired  almost  simultaneously  at  a  1.5  it  level  (or  «93  %).  Tor  comparison, 
it  can  be  noted  that  according  to  Hedlin  et  al.  (1989),  deviations  between  actual  and 
intended  delay  times  could  be  as  high  as  34%  of  the  delay  time  between  shots  for 
some  quarries  in  the  Pastern  U.S. 


Fxamples  of  Modulation  Functions 

Dahle  et  al.  (1989)  have  compiled  a  database  with  information  on  ripple-fired 
explosions  in  Fennoscandia  including  delay  times  between  individual  shots  and 
their  charges.  Figure  3  illustrates  schematically  the.  delays  and  charges  for  a 
ripple-firecl  explosion  in  this  database  (at  Dala  Kalk  AB  in  Sweden,  open  pit  neat- 
latitude  60.62N  and  longitude  15.11;).  it  consists  of  five  shots  with  25  ms  delay 
time;  the  charges  of  the  two  first  shots  are  about  one  third  of  the  three  last  shots. 

The  modulation  function  for  this  rippled  fired  explosion,  assuming  no  scatter 
in  delay  times,  is  drawn  as  a  heavy  line  in  Figure  4  together  with  results  from  simu¬ 
lations  assuming  scatter  with  a  standard  deviation  increasing  linearly  with  shot 
number  in  accordance  with  the  data  in  Figure  2.  The  average  curve  represents  a 
mean  of  1000  simulated  modulation  functions  using  a  normal  random  number  gen¬ 
erator  to  add  errors  to  the  lime  delays  of  individually  simulated  modulation  func¬ 
tions.  The  estimated  average  of  the  modulation  curve  shows  a  significantly  reduced 
second  neak  around  80  Hz  and  a  generally  increased  amplitude  level  at  higher  fre¬ 
quencies  compared  to  the  ideal  case  with  no  errors.  More  important,  perhaps,  is  the 
large  standard  deviation  of  the  average  curve.  This  effect  is  more  clearly  demon¬ 
strated  by  the  variation  in  the  10  simulated  modulation  functions  in  Figure  5.  There 
is  dramatic  variation  among  the  simulated  modulation  functions  which  also  differ 
significantly  from  the  theoretical  modulation  function  having  no  scatter  in  the  dela 
times  (heavy  line  in  Figure  4). 

The  curves  in  Figure  5  clearly  illustrate  that  significant  scatter  in  the  deiay 
times  may  also  hamper  the  identification  of  ripple-fired  explosions  using  methods 
based  on  the  detection  of  a  spectral  modulation  pattern. 


2.1.3  Number  of  Explosions  From  Autocorrelation  of  Modulation  Function 

Most  methods  to  identify  ripple-fired  explosions  that  have  recently  been  sug¬ 
gested  attempt  to  detect  modulation  or  scalloping  of  the  frequency  spectrum  of 
one  or  more  recorded  phases  (Baumgardt,  1988).  As  mentioned  earlier,  spectral 
scalloping  can  be  directly  related  to  the  time  delay  between  successive  shots  of  the 
ripple-fire,  but  does  not  provide  direct  information  on  the  number  of  shots.  Here 
we  discuss  the  possibility  of  estimating  both  delay  time  and  number  of  shots. 
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Distribution  of  Charges 


Figure  3.  Time  delays  and  charges  of  a  ripple-fired  explosion  at  Daia  Kalk  consist¬ 
ing  of  five  individual  blasts  (from  Dahle  et  a!..  1989). 
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NO  SCATTCn 


AVRIAGI: 


RELATIVE  AMPLITUDE 


SIMULATED  MOUDLATION  FUNCTIONS 


rniroiJENCY  (ii/> 


figure  5.  Examples  of  modulation  functions  for  the  ripple-fired  explosion  defined 
in  figure  3,  assuming  errors  in  scatter  of  delay  times  according  to  figure  2.  The 
nine  curves  were  simulated  using  random  number  generators  to  assign  values  to 
time  delays. 


As  shown  above,  tfie  Fourier  transform  of  the  modulation  function  of  a  ripple- 
fired  explosion  can  be  represented  as: 

N-l 

M(/w)--  ^exp(  j-uj-k-T). 


We  assume  no  or  negligible  errors  in  the  delay  times.  The  inverse  Fourier 
transform  of  the  power  spectrum  of  the  modulation  function  which  is  equal  to  the 
autocorrelation,  can  be  written  in  the  form: 

N- 1 


This  is  a  series  of  spikes  delayed  by  r  s  and  with  decreasing  amplitude  from  N  to  1. 

We  may  assume  that  the  recording,  r(l),  --of  a  phase  or  the  entire  wave  train- 
for  a  ripple-fired  explosion  is  the  modulation  function,  m(t).  convolved  with 
another  term: 

r(t)=m(t)*v(t). 


The  term  v(r)  takes  into  account  the  source  function  of  Individual  shots,  path 
effects,  and  instrument  response  by  convolution. 


We  use  a  heavily  smoothed  power  spectrum  of  the  recorded  signal,  r(t),  as  an 
estimate  of  the  power  spectrum  for  the  term  v(f).  The  smoothing  is  obtained  by  a 
standard  running  average  technique.  This  means  that: 


Ow)-  -  - -y—7—r 


is  an  estimate  (indicated  with  a  hat)  of  the  power  spectrum  of  the  modulation  func¬ 
tion.  By  taking  the  inverse  Fourier  transform  of  this  ratio,  we  get  an  estimate  of  the 
autocorrelation  of  the  modulation  function,  which  should  ideally  consist  of  a  series 
of  spikes,  delayed  by  r  and  with  linearly  decreasing  amplitudes. 


Test  on  Synthetic  Data 

licjure  6  shows  an  estimate  of  the  autocorrelation  of  the  modulation  function 
for  a  synthetic  example.  The  recording  at  the  NRDC  station  KKL  of  a  single  chemi¬ 
cal  explosion  (top  trace)  has  been  used  to  generate  a  signal  consisting  of  five  explo¬ 
sions  with  delay  time  124  ms  by  just  delaying  and  adding  five  copies  of  the  record¬ 
ing  for  the  single  explosion  (middle  trace).  This  synthetically  generated  "ripple- 
fired"  explosion  trace  has  a  somewhat  larger  P  coda  than  the  recording  for  the  sin¬ 
gle  explosion.  The  autocorrelation  for  the  modulation  function  at  the  bottom  trace 
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SYNTHETIC  EXAMPLE 


Single  chemical  explosion  at  NRKK 


l  .  I  m  i  l^i.1.4.1  *  4  ‘■l  *■* 

1IM  i*M 


!•••* 


TIMH  (S) 


figure  6.  Synthetic  example  illustrating  the  computation  of  the  autocorrelation 
function  for  the  modulation  function.  The  upper  trace  is  a  recording  for  a  single 
chemical  explosion  in  I!.  Kazakhstan  at  the  NRDC  station  NRKK  in  Karkaralinsk.  The 
middle  trace  was  obtained  by  adding  five  copies  of  the  top  trace  and  successively 
delayed  124  ms;  the  sampling  rate  of  the  data  is  250  IIz.  The  bottom  trace  is  the 
estimated  autocorrelation  of  the  modulation  function  showing  four  clear  spikes, 
with  delays  of  about  124  ms,  after  the  initial  large  amplitude. 
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shows  four  clear  spikes  after  the  zero-delayed  initial  large  spike.  Experiments 
using  this  recording  of  the  single  chemical  explosion  with  different  delays  and 
numbers  of  explosions  gave  similar  results. 


Examples  with  real  data 

The  autocorrelation  for  the  modulation  function  estimated  for  two  presumed 
ripple-fired  explosions,  one  recorded  at  the  NRDC  station  NRKK  in  Karkaralinsk  and 
the  other  at  the  high  frequency  element  of  the  NORI'.SS  array,  are  shown  in  Figure  7. 
In  both  cases,  clear  trains  of  spike-like  pulses  following  the  initial  large  amplitude 
(at  zero  delay)  can  be  seen.  The  number  of  spikes  and  estimated  delay  times  are 
indicated  in  the  figure.  At  this  point  in  time,  however,  no  independent  information 
is  available  on  the  characteristics  of  these  presumed  ripple-fired  explosions  to  ver¬ 
ify  the  estimates  on  delay  times  and  number  of  shots. 
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MODULATION  FUNCTIONS 


871 75  10:35:00.2  49.3N  72.7E  at  NRKK 
9  shots  with  30  ms  delay 


Explosion  Blaasjoen,  Norway  at  NRAO 


Figure  7.  Estimated  autocorrelations  of  the  modulation  function  for  two  presumed 
rippie-fired  explosions.  The  top  trace  is  for  recordings  at  the  NRDC  station  NRKK  at 
Karkaralinsk  and  the  bottom  is  for  recordings  at  the  high  frequency  element  of  the 
NORF.SS  array.  The  number  of  spikes  and  their  successive  delay,  indicated  in  the 
figure,  are  used  to  estimate  the  number  of  shots  and  time  intervals  between  shots. 
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3.1  SLOWNESS  ESTIMATION  WITH  INTERPOLATED  NORESS  DATA 


ABSTRACT 

In  this  note  we  explore  the  potential  of  using  the  NORESS  mini-array  to  locate 
large  teleseismic  events  with  data  resampled  to  very  high  frequencies  (1  kHz).  P- 
wave  recordings  from  15  explosions  at  the  Balapan  test  site  in  E.  Kazakh  between 
1984  and  1988  were  analyzed.  Time  delays  between  array  elements  of  the  initial  P- 
waves  are  determined  to  the  nearest  ms  by  interpolation  of  waveform  data  and 
cross  correlation.  The  time  delay  patterns  for  closely  spaced  events  show  high 
repeatability  and  correlate  to  some  degree  with  element  elevation.  The  smallest 
standard  deviations  in  azimuth  and  slowness  estimated  from  the  time  delays  (  0.26 
°  and  0.7  ms/degree  respectively)  are  obtained  only  if  the  data  from  the  outer  ring 
are  used.  The  standard  deviations  for  azimuth  and  slowness  from  f~k  analysis  of 
the  same  data  were  1.45  °  and  1.9  ins/km.  respectively.  The  average  P-wave  velocity 
below  the  NORESS  array  is  estimated  to  be  6.1  km  ±  0.7  km  by  fitting  time  delays  to 
a  three-dimensional  plane  wavefront.  Time  residuals  for  the  array  elements  from 
the  estimated  plane  range  between  j;  10  ms.  These  residuals  can  be  used  to  define 
the  shape  of  the  wavefront  of  the  initial  P-wave.  They  outline  a  smooth  surface  that 
deviates  from  a  plane  by  about  two  percent  over  one  wavelength.  The  shape  of  this 
smooth  surface  is  highly  repeatable  for  all  Balapan  events,  but  the  surfaces  for 
events  in  the  southwestern  and  northeastern  portions  of  the  test  site  appear  to  be 
tilted  about  2  °  relative  to  one  another,  which  is  about  four  times  the  theoretical 
value. 


3.1.1  Introduction 

In  this  note  we  explore  the  potential  of  the  NORESS  mini-array  for  epicenter 
locations  of  large  teleseismic  events  using  data  resampled  to  very  high  frequencies 
(1  kHz).  The  scatter  in  azimuth  and  slowness  estimates  of  large  underground 
nuclear  explosions  at  the  E.  Kazakh  testing  grounds  are  studied  using  precise  deter¬ 
minations  of  time  delays  between  array  elements.  This  approach  makes  it  possible 
to  relate  possible  errors  in  azimuth  and  slowness  estimates  to  anomalies  at  indivi¬ 
dual  array  elements. 


3.1.2  Data 

Table  1  lists  source  parameters,  reported  by  the  ISC,  of  the  explosions  in  E. 
Kazakh  for  which  recordings  at  NORESS  were  analyzed.  The  epicentral  distance 
between  the  E.  Kazakh  testing  grounds  and  NORESS  is  about  38  °.  The  events  were 
all  located  in  the  Balapan  area  and  had  magnitudes  in  the  range  mj,(/SC)= 5.5-G.2. 
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TABLE  1 


EVENT  SOURCE  PARAMETERS  (ISC) 


Date 

Time 

Lat(N) 

Long(E) 

mb 

1984 

Oct 

27 

1:50:10.6 

49.95 

78.83 

6.2 

1984 

Dec 

2 

3:19:06.3 

49.95 

79.03 

5.9 

1984 

Dec 

16 

3:55:05.1 

49.88 

78.82 

6.1 

1984 

Dec 

28 

3:50:10.7 

49.86 

78.75 

6.0 

1985 

1'eb 

10 

3:27:07.6 

49.88 

78.82 

5.9 

1985 

Apr 

25 

0:5706.5 

49.88 

78.69 

5.9 

1985 

Jun 

15 

0:57:00 .7 

49.85 

78.87 

6.1 

1985 

Jun 

30 

2:39:02.7 

49.84 

78.69 

6.0 

1985 

Jul 

20 

0:53:14.5 

49.95 

78.83 

6.0 

1987* 

Mar 

12 

1:57:17.0 

49.90 

78.83 

5.5 

1987* 

Apr 

17 

1:03:08.2 

49.84 

78.67 

6.0 

1987 

Aug 

2 

0:58:06.8 

49.88 

78.92 

5.9 

1987 

Nov 

15 

3:31:06.7 

49.84 

78.77 

6.1 

1988 

Nov 

12 

3:30:03.7 

50.08 

78.99 

6.0 

1988 

Dec 

17 

4:18:06.9 

49.89 

78.93 

6.0 

Data  for  events  marked  with  *  were  found  to  be  spurious  and  were  not  used  in  the 
analysis. 
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The  relative  locations  of  the  events  are  shown  in  Figure  1.  There  are  13  events 
southwest  of  the  fault  marked  on  the  map,  and  2  events  to  the  northeast.  Observa¬ 
tions  of  P-waveforms  at  the  UK  at  rays  (Marshall  cl  ai,  1985),  of  P-wave  amplitudes 
and  l.g  measurements  at  NORESS  (Ringdal  and  llokland,  1987),  and  of  I.g  measure¬ 
ments  at  Grafenberg  (Ringdal  and  Pyen,  1988)  suggest  that  these  two  areas  are 
characterized  by  different  geophysical  properties. 

Waveform  data  recorded  for  the  events  were  retrieved  from  the  databases  at 
the  Center  for  Seismic  Studies.  Only  data  recorded  at  the  short  period  vertical  com¬ 
ponents  (in  all  25)  were  used.  The  coordinates  (  x„yf,Z/)  of  the  elements  of  the 
NORESS  array  relative  to  the  center  element  NRAO  are  listed  in  Table  2.  The  data 
were  originally  sampled  at  25  ms  intervals,  which  corresponds  to  about  10  percent 
of  the  time  it  takes  a  plane  P  wavefront  from  an  E.  Kazakh  event  to  cross  the 
NORESS  array. 

Typical  unfiltered  traces  of  P-wave  signals  are  shown  in  Figure  2.  All  signals 
are  of  short  duration  and  have  a  very  high  signal-to-noise  ratio.  Their  amplitude 
spectra  are  also  sharply  peaked  around  2  Hz.  Because  of  these  characteristics  the 
signal  processing  for  this  study  was  based  on  unfiltered  records. 


3.1.3  Estimation  of  Time  Delays 

Data  windows  were  chosen  manually  for  each  signal.  The  beginning  of  the 
window  was  about  0.15  s  prior  to  the  pick  of  the  P-onset  at  the  center  element  of 
the  array  (NRAO),  and  the  end  of  the  window  was  set  just  after  the  first  zero  cross¬ 
ing  of  the  first  significant  peaks.  The  lengths  of  the  data  windows  varied  between 
1.6  -  2.7  s  for  the  events,  i.e.,  the  number  of  data  points,  J,  varied  between  64  and 
108. 

High  Frequency  Interpolation 

The  time  resolution  of  the  originally  recorded  data  (25  ms  between  successive 
samples)  was  improved  by  interpolation  to  1  ms  between  samples.  Standard  cubic 
splines  were  used  to  interpolate  the  data.  The  value  of  s(t)  for  a  time  t  between 
sampled  values  j-At  and  (j+l)-At  was  obtained  from: 

s(t)~A-s(jAt)\  K-s((J  \-  l)Af)  t  C-Yj+D-Yu , 


with 


A-0-i-l)-Af-f/At 

B=\-A 
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LATITUDE  (N) 

49.60  49.85  49.90  49.95  50.00  50.05  50.10 
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78.6  78.7  78.0  78.9 

LONGITUDE  (E) 


Figure  1.  Epicenters  of  the  events  at  the  Balapan  Test  Site. 


TABLE  2. 


ARRAY  CO-ORDINATES  AND  ELEMENT  CORRECTIONS 
Station  Coordinates  (km)  Correction(s) 

Code  x  y  '/.  Mean  Std 

NRAO  0.000  0.000  0.000  -0.0002  0.0015 


NRA1 

0.045 

0.142 

NRA2 

0.104 

-0.108 

NR  A3 

-0.146 

-0.034 

NRB1 

0.066 

0.318 

NRB2 

0.331 

0.025 

NRB3 

0.139 

-0.302 

NRB4 

-0.233 

-0.222 

NRB5 

-0.276 

0.158 

NRC1 

0.106 

0.683 

NRC2 

0.599 

0.337 

NRC3 

0.643 

-0.241 

NRC4 

0.204 

-0.661 

NRC5 

-0.400 

-0.571 

NRC6 

-0.692 

-0.053 

NRC7 

-0.451 

0.544 

NRD1 

0.189 

1.477 

NRD2 

1.094 

1.010 

NRD3 

1.489 

0.072 

NRD4 

1.186 

-0.905 

NRD5 

0.331 

-1.454 

NRD6 

-0.684 

-1.331 

NRD7 

-1.371 

-0.568 

NRD8 

-1.341 

0.411 

NRD9 

-0.806 

1.254 

-0.011 

-0.0008 

0.0006 

0.009 

0.0020 

0.0016 

-0.006 

-0.0034 
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ligure  2.  Examples  of  recorded  P-waves  at  some  of  the  elements  at  NORESS  from 
events  on  Aug  2  1987  and  Dec  17  1988. 


<  o  I'M*'7 

D«™  (b3  -«)-AJ2 

and  the  "second  derivatives",  Y),  derived  from  the  equations: 

At2/6  [Yj.v\AYJ.{.Yj^yYj,i~2YJ+Yj.x  j-2 . J-l 

Subroutines  by  Press  et  al.  (1986)  were  employed  for  the  computations.  The  values 
of  and  Yj  were  set  to  0.  giving  the  natural  cubic  spline,  which  has  zero  derivative 
on  both  end  points.  A  five  percent  cosine  taper  was  applied  to  the  resampled  data 
window.  An  example  of  interpolated  values  is  given  in  Figure  3. 

Waveform  Correlation 

For  each  event,  k,  we  obtained  time  delays,  r^,  of  the  P-wave  between  array 
elements  i  and  j  from  the  cross  correlation  function,  p^,  calculated  for  the  interpo¬ 
lated  and  tapered  data.  The  pIJk  was  calculated  in  a  straightforward  manner  using 
FFTs  to  increase  computational  speed. 

The  time  delay,  is  defined  by  the  time-shift,  u,  that  maximizes  the  cross 
correlation,  i.e.,  />,jk(T,jk)>Pijk(u)  for  all  u.  This  means  that  the  time  delays  were  deter¬ 
mined  to  the  nearest  ms. 

The  maximum  values  of  the  cross  correlation  functions,  i.e.,  />,jUTijk)  are  plotted 
in  Figure  4  against  distance  (see  Table  2): 

dym  ((Xf-x/Kyr-y/i 

There  is  large  scatter  in  the  data,  but  /Ujk  is  rarely  <0.95. 

Delays  were  calculated  for  all  combinations  of  array  element  pairs,  f=l . n  and 

j=l,...,n  for  each  event,  (ror  events  with  data  from  25  channels,  i.e.,  n= 25,  the 
number  of  combinations  is  300). 

Figure  5  shows  the  traces  aligned  for  two  of  the  events  using  delays  deter¬ 
mined  from  this  procedure. 


3.1.4  Analysis  of  Time  Delays 

Since  we  were  primarily  interested  in  using  the  estimated  time  delays  to  make 
precise  epicenter  determinations,  the  were  initially  compared  to  those  calculated 
from  a  theoretical  model  (P-tables  of  Jeffreys-Bullen).  We  therefore  assumed  that 
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rigure  3.  Example,  of  interpolated  data  values.  The  original  samples  are  indicated 
with  open  circles  and  the  interpolated  values  with  crosses.  The  data  points  cover 
0.125  s. 
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Figure  4.  Maximum  correlation  values  against  distance,  d,  between  element  pairs  for 
all  events  analyzed  here.  Pc  correlation  curves,  exp(  -d/1000)  from  a  linear  regres¬ 
sion.  and  exp(-d/10)  (from  Pg  data,  Mykkeltveit  el  ai.  1983)  are  drawn  for  com¬ 
parison. 
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r/fjure  5.  Traces  aligned  from  time  delays  determined  to  the  nearest  ms,  for  events 
on  Aug  2  1987  and  Dec  17  1988. 


the  observed  time  delay,  Tyk.  between  array  elements  /  and  j  for  event  k  could  be 
written  as: 


hjk  ---  r./Pk  l- Cik-Cjk  i  <  ,jk. 

The  theoretical  horizontal  slowness  vector,  p^  is  given  by: 

MPfc-sin(0k),/vcos(<>*)). 


where  the  scalar  value  of  the  slowness,  Pk,  was  calculated  from  the  travel  time 
tables,  and  the  back-azimuth,  <fk,  was  calculated  relative  to  the  center  point,  NRAO, 
of  NORESS  from  the  epicentral  coordinates  reported  by  ISC  (see  Table  1).  The  vec¬ 
tor  ?ij  is  defined  by  the  relative  coordinates  of  the  two  array  elements  (see  Table  2): 

n^(Xi-Xj,yi-yj). 


The  terms  clfc  for  array  element  i  ( /=  1,2 . n)  are  called  element  corrections  and 

are  analogous  to  station  corrections  of  travel  times  for  seismological  station  net¬ 
works.  The  error  terms  t0k  were  assumed  to  be  Gaussian  with  zero  mean  and  stan¬ 
dard  deviation  <r0k-  We  also  assumed  the  to  be  independent. 

Estimates  of  Element  Corrections 

for  each  event,  k,  one  can  estimate  the  element  corrections  c,  from  a  least 
squares  fit  to  the  equations  (in  all  301): 

C/k~ Cjk “ Tijk  ij'Pk  i'=  1  >2 24,  jmi,i+ 1  ,..,,25 

}\Cik--0 


The  average  c,k  is  given  by: 


_  1  K 


The  cik  averaged  over  five  of  the  events  (i.e.,  k- 1 . 5)  (all  located  in  southwestern 

Balapan)  for  which  data  were  available  for  all  elements  of  the  array  and  the  associ¬ 
ated  standard  deviations,  a (c/),  are  listed  in  Table  2.  The  estimated  element  correc¬ 
tions  were  quite  consistent  from  one  event  to  another.  This  is  reflected  by  the 
small  standard  deviations:  <r(cf)<l  ms  fui  14  of  the  25  elements  and  a(c,)<2  ms  for 
all  elements  except  NRD3  and  NRD8. 

The  standard  deviations  of  the  errors  were  r rok=0.7ms ,  with  virtually  no 
variation  among  events,  k.  The  contribution  in  standard  deviation  due  to  rounding 
error  (1  ms)  in  the  time  delays  is  about  0.3  ms. 
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Erroneous  Data 


Data  for  the  events  on  Mar  12  and  Apr  17.  1987,  were  not  used  to  obtain  the 
average  element  corrections  ct.  Figure  6  shows  the  residuals  (cik-ct)  for  these  two 
events  for  array  elements  i.  A  clear  offset  in  qk  of  25  ms  can  be  seen  for  several  of 
the  elements  for  the  event  on  Mar  12  and  for  one  element  for  the  event  on  Apr  17. 
The  source  of  these  apparent  sampling  errors  is  not  known,  and  the  data  for  these 
two  events  are  omitted  in  the  subsequent  analysis.  The  data  for  all  events  were 
processed  with  identical  procedures.  It  seems,  therefore,  unlikely  that  the  errors 
are  due  to  software  errors  since  sampling  errors  of  this  kind  were  detected  only  for 
2  of  the  15  events. 

Effect  of  Elevation 

The  NORF.SS  array  is  not  perfectly  horizontal.  Elevations  for  the  array  ele¬ 
ments  available  at  the  Center  for  Seismic  Studies  (cf.  /.-coordinates  in  Table  2)  vary 
by  178  m.  Assuming  a  velocity  «-G.O  km/s  in  the  surface  medium,  this  corresponds 
to  a  time  delay  of  about  25  ms  for  a  plane  P  wavefront  from  E.  Kazakh.  This  in  turn 
is  about  10%  of  the  time  it  takes  a  wavefront  to  cross  the  array  horizontally.  For 
more  distant  events  the  elfect  of  varying  surface  elevation  will  be  more  pro¬ 
nounced. 

Figure  7  shows  the  element  corrections  c,  as  a  function  of  elevation,  z,.  A  line 
(  with  the  vertical  slowness  /v- 0.183  s/km  (  pz-  |l/or?-p2j  )  is  drawn  for 

comparison.  There  is  some  correlation  in  the  data.  The  element  NRD3  at  the 
highest  point  of  the  array  has  the  largest  correction,  31  ms.  The  corrections  for 
some  other  elements  at  high  elevations  are  also  significant.  From  the  data  in  Table 
2,  it  can  be  seen  that  the  element  NRD3  (at  the  highest  point  of  the  array)  also  has 
the  largest  standard  deviation  of  its  estimated  element  correction,  a(ct),  among  the 
elements  (see  Table  2). 

Factors  other  than  elevation  clearly  influence  the  element  corrections,  as  indi¬ 
cated  by  the  five  outlying  data  points  in  Figure  7.  In  Figure  8  the  corrections  are 
plotted  as  a  function  of  relative  horizontal  location  of  the  elements.  A  correction 
for  delay  due  to  elevation,  z,-pz  with  the  vertical  slowness  pz=0.183  s/km,  has  been 
applied  to  the  c(,  -  values.  Positive  and  negative  corrections  are  separated  by  a  line 
running  in  the  northwest-southeast  direction.  There  is  no  apparent  correlation  with 
the  local  surface  geology,  some  features' of  which  are  also  indicated  in  the  figure. 

Scatter  in  Azimuth  and  Slowness  Estimates 

We  have  looked  at  time  delays  among  elements  relative  to  values  calculated 
from  theoretical  slowness  and  azimuth  values.  The  reverse  problem  of  estimating 
slowness  and  azimuth  is  of  more  practical  interest.  Such  estimates  can  be  obtained 
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Figure  8.  Average  element  corrections  (with  correction  for  elevation)  as  a  function 
of  horizontal  element  coordinates.  Positive  and  negative  corrections  are  indicated 
with  circles  and  squares  and  imply  "slow"  and  "fast"  delays,  respectively.  The  larger 
the  symbol  the  larger  its  absolute  value. 
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ill  a  straightforward  manner  by  fitting  the?  time  delays  to  those  of  a  plane  wave* 
front.  Because  of  the  significant  values  of  (lie  element  corrections,  C/(,  such  esti¬ 
mates  will,  in  all  likelihood,  lie  biased.  However,  by  subtracting  the  estimated  ele¬ 
ment  corrections  c,  from  the  observed  time  delays  r iJt,  unbiased  estimates  may  be 
obtained.  1'hat  is  to  say,  we  can  estimate  the  components  of  the  slowness  vector, 
Ih-  (Pitk'Pyk)  for  event  k  from  a  least  squares  lit  to  t he  equations: 

i'ijk  Pk  ■•••  Tiju  (( i  t))  i“  1  ,...n- 1 ,  jmi+ 1 ,...,n 

I  stimates  of  azimuth.  </»*..  and  slowness,  /\,  are  obtained  in  t lie  usual  manner: 

ota»{)\k/Pyk)< 

Pk  (pu  i  Kkf- 


f.rrors  in  azimuth  and  slowness  are  defined  as  the  difference  between  estimated 
values  and  those  calculated  from  the  Jeffroys-Bullcn  travel  lime  model  and  the 
coordinates  reported  by  the  ISC:  * 

A>!>k  <Ac-  ‘Ac. 


^Pk  -  P c  Pk- 


The  quality  of  such  unbiased  estimates  may  then  be  described  by  standard 
deviations  of  the  errors,  and  a(A p),  for  t lie  events.  The  element  corrections 

estimated  from  data  of  five  of  the  events  in  the  southwestern  part  of  the  Balapan 
test  site  are  applied  to  the  time  delays  for  all  events.  As  the  estimates  of  element 
corrections  for  NRD3  and  NRD8  have  large  standard  deviations,  time  delays  involv¬ 
ing  those  elements  were  omitted  in  the  estimation. 

We  estimated  the  slowness  and  azimuth  for  each  ring  of  the  array  separately. 
I  stimates  from,  for  example,  the  B-ring  were  thus  based  on  time  delays  among  the 
elements  NRB1,  NRB2,  NRB3,  NRB4,  and  NRB5.  The  II  events  in  the  southwestern 
part  of  Balapan  were  analyzed  first,  figures  0  and  10  show  the  standard  deviation 
of  t he  azimuth  and  slowness  errors,  respectively,  for  the  four  rings.  The  standard 
deviation  in  azimuth  is  reduced  from  just  below  2.0*  for  the  A-ring  to  0.20*  for  the 
D-ring.  The  standard  deviation  in  the  theoretical  azimuth  for  these  1 1  events  is 
0.0K“.  The  standard  deviation  in  slowness  is  reduced  in  a  similar  way  from  approxi¬ 
mately  0.004  s/km  to  0.0007  s/km.  The  smallest  standard  deviation  in  slowness 
corresponds  to  a  standard  deviation  in  estimated  epicentral  distance  of  about  1.45*. 
It  can  be  noted  that  the  standard  deviations  for  the  entire  array  (except  elements 
NRD3  and  NRD8)  are  somewhat  larger  than  those  for  the  D-rlng.  This  means  that 
tin*  best  results  lor  estimating  slowness  and  azimuth  of  this  group  oT  events  Is 
obtained  from  relative  time  delays  of  7  of  the  25  elements  of  the  array. 
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figure  9.  Standard  deviation  of  azimuth  error  for  1 1  events  in  southwestern 
Balapan  as  a  function  of  array  tint',. 


3- 


STD  (Slowness)  s/km 


0.005 


Aperture  (km) 


/  iinuc  10.  Standard  deviation  of  slowness  error  for  II  events  in  southwestern 
Balapan  as  a  function  of  array  ring. 
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f  igure  1 1  shows  the  azimuth  error  as  a  function  of  theoretical  azimuth  for  the 
D-ring  estimates.  We  have  included  the  data  for  the  two  events  in  the  northeastern 
Balapan.  These  two  data  points  fall  in  the  lower  left  corner  of  the  diagram,  which 
corresponds  to  a  shift  of  1.7“  for  a  change  in  theoretical  azimuth  of  about  only  0.3*. 
This  shift  in  azimuth  is  caused  by  shifts  in  the  element  corrections,  C/*,  for  the  two 
northeastern  Balapan  events  relative  the  the  average  corrections,  c,,  for  the 
southwestern  Balapan  events,  figure  12  shows  the  residuals,  c^-C/,  for  the  two 
events.  Positive  and  negative  shills  of  about  ms  can  he  noticed. 

Comparison  with  f-k  analysis 

For  the  sake  of  comparison,  azimuth  and  slowness  values  were  also  calculated 
using  f-k  analysis,  which  has  become  a  standard  processing  techniques  for  small 
aperture  seismic  arrays.  The  algorithm  based  on  broad  band  f-k  analysis  at  the 
Center  for  Seismic  Studies  was  employed  for  this  purpose.  Data  for  all  elements 
were  used  (except  for  channels  witli  spurious  data).  The  f-k  analysis  was  applied 
to  the  same  data  windows  as  used  in  the  correlation  analysis  here.  The  standard 
deviations  in  the  azimuth  and  slowness  errors  for  the  11  events  in  the 
southwestern  part  of  the  Balapan  test  site  were.  1.45°  and  0.0019  s/km,  respectively. 
This  azimuth  error  is  comparable  to  that  of  the  B  ring  data  and  about  five  times  the 
D-ring  error  in  figure  9. 

Planarity  of  the  wavefront 

Up  to  this  point  we  have  compared  the  time  delays  with  values  corresponding 
to  a  plane  wavefront  derived  from  a  theoretical  velocity  model.  The  time  delays  can, 
of  course,  also  be  used  to  estimate  the  slowness  and  azimuth  of  a  plane  wavefront 
without  applying  any  element  corrections.  Such  estimates  will  in  all  likelihood  be 
biased,  but  it  is  of  some  interest  to  see  how  well  the  time  delays  fit  those  of  an 
estimated  plane  wavefront. 

Because  of  the  differences  in  elevation  among  the  array  elements,  we  estimate 
the  components  of  a  three  dimensional  slowness  vector  (f>k-(Pxk>Pyk>Ptk))  for  event 
k  from  a  least  squares  fit  to  the  lime  delays: 

fi,  Pk  rijk  0“  1 . n-1;  j-l+1 . n) 

ITrors  in  the  estimated  values  for  azimuth  and  horizontal  slowness,  i.e., 
atan(px/py )  and  .  for  the  I  I  events  at  southwestern  Balapan  were  0.28  * 

0.80“  and  -0.0088  ±  0.0014  s/km,  respectively.  The  vertical  slowness  pt  is  related 
to  the  surface  P-wave  velocity,  a,  and  the  horizontal  slowness  by: 

-v  P*  I  Py  t  P/- 
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figure  12.  Difference  between  element  < 
Balapan  (Dec  2  1984  and  Nov  12  1988)  ar 
southwestern  Balapan.  The  differences  are 
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Estimated  mean  and  standard  deviations  of  «»  from  tills  relation  for  the  11 
southwestern  Ralapan  events  were  6.1  i  0.7  km/s. 

The  residuals.  cik,  from  the  estimated  plane  wavefront  at  element  /  for  event  k 
can  be  estimated  in  a  way  similar  to  the  element  corrections  above: 

Cik-Cj^-rijk  tljrk  (l-l . n-1; j-i+1 . n) 


here  pk  represent  tlie  3-dimensional  slowness  vector  already  estimated  from  the 
time  delays. 

The  average  residual  for  the  1 1  events  in  southwestern  llalapan  has  a  standard 
deviation  of  u(c/)-4.5m.v.  l  or  comparison,  the  standard  deviation  for  data  estimated 
by  fitting  a  plane  wave  with  only  a  horizontal  slowness,  i.e.,  only  px  and  py  com¬ 
ponents  of  p,  was  6.3  ms. 

Average  residuals  for  the  I  I  events  in  the  southwestern  part  of  Ralapan,  which 
range  between  -10  to  10  ms,  are  shown  as  a  function  of  horizontal  location  in  Fig¬ 
ure  13  in  a  perspective  diagram.  If  the  velocity  below  NORESS  is  about  6  km/s,  this 
corresponds  to  an  undulation  of  about  60  in  of  the  wavefront  over  a  horizontal  dis¬ 
tance  of  about  2.67  km  (about  one  wavelength),  which  corresponds  to  about  2%. 
Positive  and  negative  residuals  divide  the  plane  into  two  portions:  the  central 
northern  and  southern  part  of  the  array  have  positive  residuals,  whereas  the  east 
west  and  central  parts  have  negative  residuals.  The  residuals  define  the  surface  of 
the  front  of  the  initial  P-wavc. 

The  residuals  are  also  plotted  in  the  horizontal  plane  in  Figure  14,  where  some 
features  of  the  surface  geology  at  the  array  are  indicated  (Mykkeltveit,  1987).  Some 
of  the  "slow"  portion  of  the  wavefront  coincides  with  a  large  tongue  of  gabbro 
extending  almost  due  east  from  the  center  of  the  array. 

The  residuals  for  the  two  events  in  the  northeastern  pari  of  Ralapan  are  very 
similar  to  those  of  the  southwestern  part.  However,  the  errors  in  azimuth  and  slow¬ 
ness  arc  clearly  different,  *1.73  "  and  -0.0066  s/km,  respectively.  Figure  15  shows 
the  direction  of  the  normal  to  the  estimated  plane  wavefronts  (assuming  the  surface 
velocity  <i-6  km/s)  projected  on  the  lower  hemisphere.  It  appears,  therefore,  as  if 
the  shape  of  the  wavefronts  from  the  two  areas  of  Ralapan  are  very  similar,  but  they 
are  slightly  tilted  relative  to  one  another  differing  by  1.55*  In  backazimuth  and 
0.82°  relative  to  the  vertical  (incidence  angle).  * 

It  might  be  more  plausible  to  relate  this  tilt  to  the  local  structure  beneath  the 
NOItrSS  array  than  to  efiects  near  the  source  regions  or  along  the  propagation 
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ritjinc  14.  Plane  wave  residuals  as  a  function  of  horizontal  coordinates.  Circles  and 
squares  are  used  for  positive  and  negative  values,  respectively,  corresponding  to 
"slow"  and  "fast"  portions  of  the  wavefront.  Some  surface  geological  features  (after 
Mykkeltveil,  1987)  are  indkatcd. 
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figure  15.  Lower  hemisphere  project  ion  oi  the  directions  of  a  theoretical  plane 
wave  and  apparent  plane  waves  for  the  two  portions  of  the  Balapan  test  site. 


paths.  However,  l here  are  observations  from  teieseismic  recordings  at  other 
seismic  stations  dial  show  dear  differences  for  the  northeast  and  southwest  sec* 
lions  of  Bainpan.  Systematic  differences  in  short  period  waveforms  at  the  UK  arrays 
were  noted  by  Marshall  cl  al.  (1985),  who  also  suggested  that  these  differences 
might  be  due  to  systematic  differences  in  corner  frequencies  of  explosions  at  the 
two  regions.  Analyzing  broad  band  recordings,  Stewart  (1988),  however,  suggests 
that  the  differences  are  more  likely  to  be  related  to  source  depth  of  the  explosions. 
Ringdal  and  Marshall  (1989)  suggested  that  the  two  sections  of  Balapan  are  charac¬ 
terized  by  different  geophysical  properties.  They  observed  a  systematic  variation  in 
Him  I'  I  n  I »l*m  li»i‘  mkpIumiimh  nn-imt  Hit*  Ihilripmi  ii'hiihi.  Fuf  ttvtMiUt  with  similtti'  U 
magnitudes,  the  P  amplitudes  from  explosions  in  the  northeast  were  smaller  than 
those  from  events  in  the  southwest.  If  the  tilt  observed  here  is  attributed  to  the 
near  source  structure  it  appears  that  the  outgoing  rays  are  deflected  more  strongly 
from  events  in  the  northeast  than  in  the  southwest  part. 


3.1.5  Concluding  Remarks 

In  this  study,  we  analyzed  NORF.SS  P-wave  recordings  for  15  explosions  from 
the  Balapan  lest  site  in  F.  Kazakli  between  1984  and  1988.  P-wave  time  delays 
between  array  element  pairs  were  determined  to  the  nearest  ms  from  cross  correla¬ 
tions  of  interpolated  waveform  data.  The  pattern  of  time  delays  shows  high  repea¬ 
tability  for  events  in  the  same  area  of  Balapan. 

Kslimates  of  time  delay  corrections  for  the  elements  relative  to  delays  for  a 
theoretical  plane  wavefront  varied  between  -19  to  31  ms.  The  standard  error  in  the 
estimation  was  0.7  ms.  The  element  corrections  correlate  to  some  degree  with  ele¬ 
ment  elevation. 

Azimuth  and  slowness  were  estimated  from  the  lime  delays  and  the  element 
corrections.  The  best  results  were  obtained  if  data  from  only  the  outer  ring  of 
NOR RS. S  were  used.  The  standard  deviations  in  azimuth  and  slowness  estimates, 
0 .21)  "  and  0.7  ms/degree.  respectively,  in  this  case  were  significantly  smaller  than 
those  obtained  from  f  -  A:  analysis  on  the  same  data,  1.45  *  and  1.9  ms/km,  respec¬ 
tively. 

The  average  P-wave  velocity  below  the  NORF.SS  array  was  estimated  to  be  6.1  ± 
0.7  km/s  by  filling  the  time  delays  to  those  of  a  three-dimensional  plane  wavefront. 

Flement  travel-time  residuals  from  the  estimated  plane  wavefront  varied 
between  ±  10  ms.  These  residuals  define  a  smooth  surface  for  the  wavefront  of  the 
initial  P-wave.  The  maximum  undulation  of  this  surface  was  about  60  m  or  about  2% 
ac  ross  one  wavelength  of  the  plane.  The  shape  of  wavefront  surface  was  highly 
repeatable  for  all  Balapan  events,  whereas  the  wavefronts  for  events  in  the 
southwestern  and  northeastern  Balapan  appeared  to  be  tilted  about  2  *  relative  to 
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one  nnolher,  which  is'  about  four  times  the  theoretical  value.  This  implies  that  out¬ 
going  rays  from  the  Balapan  area  may  be  deflected  more  strongly  in  the  near  source 
region  from  events  in  the  northeast  section  of  the  test  site  than  from  events  in  the 
southwestern  section. 


Hans  Israelsson 
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Los  Alamos  National  Laboratory 

P.O.  Box  1663 
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Los  Alamos,  NM  87545 
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DARPA/NMRO 

1400  Wilson  Boulevard 

Arlington,  VA  22209-2308 
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Los  Alamos  Scientific  Laboratory 

P.  O.  Box  1663 

Mail  Stop  D-406 

Los  Alamos,  NM  87545 

Dr.  Michael  Shore 

Defense  Nuclear  Agency/SPSS 

6801  Telegraph  Road 

Alexandria,  VA  22310 
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United  States  Senate 

Washington,  DC  20510 
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Dr.  Larry  Turnbull 
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Cameron  Station 
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Washington,  DC  20330 
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Arlington,  VA  22209 
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Hanscom  AFB,  MA  01731-5000 
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Hanscom  AFB,  MA  01731-5000 
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1400  Wilson  Boulevard 

Arlington,  VA  22209 
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